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Frequency integrated radiation models for 
absorbing and scattering media 

By J.F. Ripoll and A.A. Wrayf 


1. Motivation and objectives 

The objective of this work is to contribute to the simplification of existing radiation 
models used in complex emitting, absorbing, scattering media. The application in view 
is the the computation of flows occuring in such complex media, such as certain stellar 
interiors or combusting gases. In these problems, especially when scattering is present, 
the complexity of the radiative transfer leads to a high numerical cost, which is often 
avoided by simply neglecting it. This complexity lies partly in the strong dependence of 
the spectral coefficients on frequency (Modest 2003; Siegel & Howell 2001). Models are 
then needed to capture the effects of the radiation when one cannot afford to directly 
solve for it. In this work, the frequency dependence will be modeled and integrated 
out in order to retain only the average effects. A frequency-integrated radiative transfer 
equation (RTE) will be derived. In it, the absorption and scattering will be treated 
through the use of mean coefficients (Siegel & Howell 2001 and references in it). To 
obtain these coefficients, it is needed to assume a form for the intensity, which we take 
to be the maximum entropy closure (Minerbo 1978). Such an intensity is a function of 
the macroscopic radiative energy and flux and accounts for the variations of radiation 
in the considered medium. Models for mean absorption, mean isotropic scattering, and 
mean non-isotropic scattering coefficients will be proposed in the case where the various 
spectral coefficients can be written as polynomial functions of the frequency. Some of 
these models have already been derived and tested for non-scattering media in (Ripoll 
et al. 2001; Ripoll & Wray 2004a). They are here extended to the general case of emitting, 
absorbing, and scattering media. A direct application will be given for soot, which follows 
a linear frequency law for absorption; isotropic and incoming scattering spectral models 
are also here roughly approximated by a linear law. Macroscopic radiation models will also 
be derived with absorption and scattering coefficients since they constitute an alternative 
to the use of the RTE in cases where this equation is too costly to solve, as in many 
coupled problems. Finally, we believe another application of these models could be for 
radiating flows occuring in dusty media. 


2. A frequency integrated RTE with mean coefficients 

2.1. Generalities 

The radiative transfer equation (RTE) describes the evolution of the radiative intensity 
within a emitting, absorbing, and scattering medium and is given by 

-dtl+n-VI = a a MB(iy,T)-a a (iy)I-a is Ml+^^ [ ^ -> n')dn' (2.1) 

c 4 tt J n , 
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where the intensity I = I(t, r, SI, v) is a function of the time t, the position r, the direction 
of propagation SI, and the frequency v. Here c is the velocity of light, o a (v) is the 
spectral absorption coefficient, a ls (v) the spectral isotropic scattering coefficient, cr s (v) 
the spectral incoming scattering coefficient, and O' the original direction of radiation 
scattered into SI. The Planck radiative intensity B describes the isotropic emission of the 
medium at the frequency v and temperature T by 

B (,,r) = ^[e x p(^)-ir‘ (2.2) 

where h is the Planck constant, k the Boltzmann constant, and v the frequency. We 
assume in this paper that the phase function of scattering $(v, SI — > SI') can be expressed 
as 

$1 — > O') = a{u) + ffiv)Sl ■ St' + 7 (v)(fi • St') 2 + r](v)6(Sl - Si') (2.3) 

where S(-) is the Dirac delta function, and the coefficients, a, B. 7 , and r] are all in general 
functions of the frequency. They must be defined such that the following normalization 
property holds: 

f $(v,Si^Si')dSt = a is (v). (2.4) 


2.2. Derivation of the RTE with mean coefficients 

By integrating the RTE over frequency and introducing the quantity! J(t, r, St) 
= J 0 °° / (t, r,Sl,v)dv = {/)„, we obtain 


- c d t J + SI ■ VJ = (cr a {v)B{y, T))„ - {cr a {v)I) v - {a s (v)I)„ 

+ ( ~T~^ f /(n')$(v, ^ Si')dSt' )„. 

^ Jn< 

We now introduce the following mean absorption and scattering coefficients: 


'CO 


{o a (v)B(v,T)) 

(B(v,T)) v 


(2.5) 


(2.6) 


<7“- v (t,r,n) = 


{a a {v)I{t,T,Sl,v)) v 


a“-"(t,r,n) = 


E (I(t, T,Sk,u)) v 

- a s {v)r)(v)/ 4tt] 7(f,r, St,v)) v 


{ I(t,T,Sl,l/)) v 


_ s u 

a E\t, r ) ~ 


a|r(t,r, SI) = 


(I(t, r, Sl',v))„w 
{a s (p)f3(o)Sl ■ Si'Ijt, r, Sl',v)) VtSV 
{St-St'I(t, r, n',i/)) Vi n- 


(2.7) 

(2.8) 

(2.9) 

(2.10) 


0c(Cr, «) 


(o- s (i/) 7 (^)(r? • Si') 2 I(t, r, nv)) y; n- 
((St-Si') 2 I(t, r, fi', «/)>„,!,< 


( 2 . 11 ) 


where the superscripts e, a, is, and s designate respectively emission, absorption, isotropic 
scattering, and scattering. We define x{ v ) = cr s (v)r](i') / (4na ts (v)) . Using these defini- 
tions and without any assumptions other than the form of the phase function (2.3), the 


t We will denote the integration of a function / over the variables X , Y . Z as (f)x,Y,z 
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frequency-integrated RTE becomes 


1 


dt J + • V J — 


-f—aT 4 - 

4n 


'J + 




+a s F n-n l +(r s G (n-n 1 ) 2 ) dn' 


(2.12) 

where the constant a = 8n 5 k 4 /(15h 3 c 3 ). It should be noted that the mean scattering 
coefficients in the integral do not depend on f l, unlike ofj v and a l E ,v . 

As a simplifying approximation, we now eliminate the O dependence of the two mean 
coefficients a E and a F . This can be done by simply approximating the numerator and 
denominator in (2.7) and (2.8) by their ^-integrated forms. Alternately, we integrate 
2.12 indefinitely over f l , using the polar variables n = cos 9 and <f > , to obtain 


n 


a , j + s! • vjdtd* = ?faT^ - 

C 47T 


mt, r,n, v)) v dfid<t> 

f{a is (v) [1 - xW] I(t, r, J2, v)) v d/j,d<f> 


/ 


/ 


J dfid(f> 
J d/idcf) 


+ 


f(I(t,r,S2j v)) v dgi d<p 

[ - 1 - [ j(n')(a s E + a s F n-n' + a s G (n-n') 2 )dn' d/j,d<t> 

J 4tt Jfii 

(2.13) 


We now replace the indefinite integrals in the numerator and denominator of the second 
and third terms on the rhs of (2.13) with definite integrals over the full 4n of fl and then 
differentiate the resulting equation with respect to jj, and <f>. 

This approximation allows us to define the following new mean absorption and scat- 
tering coefficients a F , a E , and a E , the last two being approximations of a E v and a % E v 
respectively: 


Vp(T) 

o a E (t, r) 

vstf, r) 


{o a (v)B(v,T)) v & _ e ,„ 

{B(v,T)) Vttl ap ( j 

{a a (v)I(t, r, n, v)) v ,n _ {a a (v)I(t } r,^l,v)) uM 

{lit, r,il,v)) v ,n ~ {I(t,r,{l,v))„'4* 

{a is {v) [1 - xM] v)) v ,n 

{I{t,T,Sl,v)) Vt n 


(a is (v) [1 - xi v )\ v)) v ,i tt 


(2.14) 

(2.15) 


(2.16) 


The elimination of the il dependence of the mean absorption coefficients constitutes our 
first approximation. It is not strictly necessary to the closure: one could choose to not 
make it. The three mean incoming scattering coefficients are unchanged from (2.9)-(2.11). 
The frequency-integrated RTE with these approximations becomes 


-dtJ + ^ • VJ = ^-aT 4 - a%J - cr E J 
c 47 r 

+ - 1 - / J(n')(a' E + a s F n- n' + a s a (n- n') 2 )dn' . (2.17) 

4 ?r Jn< 


We now introduce the following macroscopic quantities: the radiative energy given by 


E R (t, r) 


;W, r,n,i/))„ jn = ^{J{t,rSl)) a 


(2.18) 
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F R (t,r) = r,fi)) n (2.19) 

and finally the radiative pressure 

P n(t, r) = ^ (f2 ® fll(t,r, fi, v)) v ^ = ^ (n <g> flJ(t,r, n)) n (2.20) 

It should be noted that in (2.17) these moments appear in the integral term of the 
scattering. This term can then be absorbed in the definition of the moments and the 
integrated RTE written as 

-dtJ + n ■ VJ = ^«T 4 - a%J - a'WJ 

C 47T 

+ ^(ca s E E R + a s F F R -n + ca s G ^ PgSlW). (2.21) 

n i,j= 1-.3 

For this equation four remarks should be given. First, the introduction of the mean coef- 
ficients allows expressing the scattering term as a function of the moments. Integrations 
of the intensity are still needed in order to compute the moments, and in that sense the 
formulations (2.17) and ( 2 . 21 ) are equivalent. Second, if the phase function has moments 
of order higher than 2 , this will introduce moments of order higher than the pressure in 
the scattering term. Third, if the intensity used in the mean incoming scattering coeffi- 
cients is the exact one, then the incoming scattering term is exact. Finally, this equation 
has been derived with two assumptions relative to the absorption and isotropic scattering 
terms, namely (2.15) and (2.16). 

The frequency-integrated RTE (2.21) will be closed in the next section by propos- 
ing expressions for the mean coefficients Op, a%, a' F . a%, a s F , and a G in terms of the 
microscopic spectral coefficients a a , a ts , a s , a, /?, 7 , and r]. 


2.3. Closure of the radiative equation with mean coefficients 
The mean coefficients (2.15), (2.16), (2.9), (2.10), and (2.11) are closed in this section by 
assuming a particular functional form for the intensity used in their definitions. For a% 
for instance, it is assumed that 


a = (a a (^)/(t, r,n, */))„, 47T _ {a a (n)I*(t, r, », i/))„, 47r = a * 
° E {i{t, r,n,v))„ i4n ~ {i*{t,r,n,v))„ i4ir aE 

where the assumed intensity is denoted I*. Similarly we assume that 

As is* S _s * s S * S _s * 

&E — a E 5 a E — a E > a F — a F , a G _ a G 


(2.22) 


(2.23) 


where the designation * for the mean coefficients, o'f* . a F *, a%*, a F *, and a G *, indicates 
that the intensity I has been replaced by the pseudo-intensity I*. 

If we now assume that, at a microscopic level, the spectral absorption coefficient o a can 
be approximated as a sum of polynomial functions of frequency: o a {v) = Cft' 1-3 , 

it follows that 


JV N 

a c p=Y, a A c ?-T) and a%* = £ a* Ei (C?, T R , f) (2.24) 

*= 1 i = 1 


where ap i and OEi will be derived in the next section (Eqs. (2.32) and (2.38)). In par- 
ticular, it will be shown that they depend on the radiative temperature T R and on the 
anisotropic factor f defined below. 
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If the spectral scattering coefficient can also be written as a zs (n) [1 — x(4 = 

Sill C-V*~ 3 , then the isotropic mean scattering coefficient can be then written in terms 
of UEi (see Eq. (2.38) in the next section) and is given by 


N 

4 s * = E CT S(C'f,T R ,f) (2.25) 

i = 1 

Now let us assume that the incoming spectral scattering coefficients can also be approx- 
imated as sums of polynomial functions: 


N 


N 


N 


a s (n)a(n) = ECT’V -3 . a s (u)/3(u) = Ecf-V” 3 , and a s {v)l{v) = E^4 


i 3 


i~l 


i = 1 


i = 1 


(2.26) 


then we have 

4 * = jt<r*Ei(C?',TR, f), 4 * = E4 i (C'f’ S 4i J , f), 4 * = E 4i(^f S 4, f). 


N 


N 


i—1 


i—1 


i—1 


(2.27) 

Finally, using these models for the mean absorption coefficients, the frequency-integrated 
RTE in its closed form is given by 


-dtJ + • VJ = ^«T 4 - a%*J - a%*J + 

c 47 r 


l-(ca s E *E R +a s F *F R -n + ca s G * E P rWW) (2.28) 


ij= 1-3 


where all the mean coefficients are defined in terms of the functions a% v a* Fi , and a q 
defined in the following section. 


2.4. Computation of the mean coefficients 


The pseudo-intensity which is used in the definition of the mean coefficients is obtained 
from the maximization of the radiative entropy (Minerbo 1978; Fort 1997) and is given 
by 


I*(t, r,n,i/) = 


2 hn s 


exp( 


hv 


kT*(n) 


)-l 


with T*(Sl) = — A • fi)) and A and B defined by 


(2.29) 


_ 2-v/4-3||f|P 

||f|| 2 


B T r 


3+IIAH 2 

3(1 — 1 1 A || 2 ) 3 


(2.30) 


These two coefficients are defined from the macroscopic quantities Tr, the radiative tem- 
perature, and f, the anisotropic factor. The radiative temperature is defined in terms of 
the radiative energy by Er = aT and the anisotropic factor is given by f = F r/(cEr). 

Let us assume that a polynomial approximation in frequency of the spectral absorption 
and scattering coefficients can be done. Such a coefficient is chosen to have the following 
form 

N 

a (^) = E C i l ' t 3 

i = 1 


(2.31) 



286 


J.F. Ripoll AND A. A. Wray 

where C, are constants which depend on the gaseous medium considered, the volume 
fraction of the main species, the pressure, etc. 

First, the well-known Planck mean absorption coefficient is computed from the Planck 
function (see Ripoll et al. 2001): 


ap 


(T) = 


(B( V ,T)) v>a 


N 


15 


N 


= J2°n(Ci,T) = - ^i\ Ci ai + l) 


i—1 


i—1 


kT 


i - 3 


Second, the three mean coefficients &E; and gq are defined by 

{a{v)r(t, r, 12,i/))„, n 


a* F {t,r, 12 ) = 


{g(p)£1 • 1 l'I*(t, r, 12', v)) v <n> _ 12 • (g(v)£1' I*(t, r, 12', v)) v> n' 


(12 • 12 'I*(t, r, 12', 12 • (12'/*(t, r, 12', v))u,ri' 

Gait, r, 12) = r,12',i/)) v , n . 


(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 


((12 • 12') 2 J*(f,r, 12', v)) v ^w ' 

It should be noticed here that in I* the 12 dependence is only present through a 
scalar product with the vector A. This greatly simplifies gf as follows. The vector 
(12'/*(t, r, 12 ' ,v))si' necessarily has the form X(t, r, v) A (t, r), where X is a scalar, which 
leads to 

{a{v)X{t,r,v)) v 


G* F {t, r) = 


(2.37) 


{X(t, r, v)) v 

Thus, Gp does not depend on 12. 

Moreover, the simple form of I* leads to analytical expressions for ge and g Fi though 
the computation is not detailed here (see Ripoll et al. 2001). These coefficients have also 
been tested and validated for simple ID problems (see Ripoll & Wray 2004a). They are 
given by: 

N 


a* E (T R ,f) = ^a* Ei (Ci,T R ,f) 


i—1 


45 


N 


k 


hB( 1 - || A|| 2 ) 


= 2^ £<* - *)! C. C(« + 1) 

i = 1 
N 

G* F (T R ,f) = J2°h(Ci,T R ,f) 


i - 3 


^(l|A| 


(2.38) 


i—1 


45 

8n 4 


i - 3 


^(1|A| 

II All 3 


with 


Pi(||A||) = ((1 + ||A||) 4 - (1 - || AH)*)/ (3 + || A|| 2 ) 
P/(||A||) = (1 - ||A||)*(i||A|| + !) + (! + ||A||)‘(*||A|| - 1) 


(2.39) 


(2.40) 

(2.41) 


with ( the Riemann Zeta function and A and B given in (2.30)f. The derivation of gq 


f for i real but non-integer, these expressions are valid provided i\ is replaced with T(« + 1) 
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involves the pressure tensor and is much more complex. We have 


°G = 


j,k = 1..3 


(52j,k £ 2 ' 

Following the same developments as in Ripoll et al. 2001, it is found that 

N 

a* G (T R ,f) = J2^G i iC i ,T R ,f) 


i = 1 


4c 

= 5TI X>- 3 > !c, .« i + 1 ) 


N 


2tt 4 


i = 1 


k 


i - 3 


^(a||A| 


where P G ( A, ||f2||) is given by 

pi m 

G 1 - ||A|| 2 + 4(A - O) 2 


(2.42) 

(2.43) 

(2.44) 


+ 


(1 + IIAID* 


1 — || A|| 2 + 4(A • siy 


hB{\ - ||A|| 2 ) y II A|| 6 

[(1 _ II A|| 2 ) || A|| 2 + (|| A|| 2 - 3)(A • O) 2 
+*I|A|| [(1 + ||A||)||A|| 2 - (|| A ||* + 3)(A • f2) 2 ]] 
■[(||A|| 2 - 1)||A|| 2 + (3 - ||A|| 2 )(A - O) 2 ) 

+i||A|| [(1 - ||A||) || A|| 2 + (||A||t - 3) (A • O) 2 ]] (2.45) 


It should be noticed here that the Planck mean a R is a function of the volume fractions 
Ci and of the temperature, while the effective mean coefficients a* E . a F are functions 
of Ci and of both the radiative temperature and anisotropic factor. In addition, a G is 
dependent on the scalar product between a vector parallel to the flux and the direction 

f2. 

The mean coefficients have been introduced to avoid the cost of the frequency inte- 
gration. Thus, to be useful, the numerical cost of evaluating the mean coefficients and of 
the iterations needed to solve for J must be lower than that of these time integrations. 
The coefficients which are proposed here are analytic and should have low computational 
cost. Furthermore, the three coefficients cr EFG have common parts (see (2.38)-(2.44)). 
For absorbing media, it has been found that fewer than 6 iterations to solve the RTE were 
needed in the cases studied (Ripoll & Wray 2004a), which indicates that this method 
does not have a high numerical cost. 


3. The linear case: application to soot 

The apparent complexity of the coefficients in the previous section is only due to the 
polynomial expansion in frequency. For sooty media the absorption dependence is linear, 
and in this case the coefficients turn out to be quite simple. 

The linear spectral absorption and isotropic scattering coefficients are written as 

o e {v) = C e v , a a {v)=C a v , and a is {v) [1 - xW] = C is v (3.1) 

As a rough approximation, we here also assume that the scattering follows a linear lawf. 
The spectral scattering coefficients are then written as 

a s (n)a(n) =C a ’ s v, a s (n)/3(n) = C p ’ s v, and a s (n)j(n) = C 7 ’V. (3.2) 

f More complex models, such as in Houf 1999, could also be treated but would need to 
introduce higher order terms in frequency 
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We define for convenience the following coefficients: 

Cp = 360 ^-C x 
F n 4 h 


(3.3) 


where x is a variable which will alternatively take the value a, e, is, (a, s ), (/3, s ), (7, s) 
below. The Planck absorption coefficient becomes then 


a e P (T) = a P4 (C e , T) = 360 C e ^-T = C e P T 

7 r 4 /i 


(3.4) 


Using (2.38), (2.39), (2.44) for i = 4, and the definition (3.1), the general form for the 
mean coefficients is given by 

f) = 3 C*p — =C* P T R G E (f ) (3.5) 


a* F4 {C x ,T R ,i) = 


2?(1 - ||A|| 2 )(3 + ||A|| 2 ) 

= C X P T R G F ( f) 


C X P 5 + 1 1 A || 2 


a* G4 {C x ,T R ,f) = C X P 


4 B{\ — || A|| 2 ) 

(1 - ||A|| 2 +6(A - O) 2 ) 


(3.6) 

= C x P T R G G (f,U) (3.7) 


^(1 — II A || 2 ) (1 — || A || 2 + 4(A ■ J2) 2 ) 

where the functions Ge,f,g are obtained by using the definitions (2.30) in (3.5)-(3.7). 
Hence, the mean absorption and isotropic scattering coefficients are given in the linear 
case by 


<J%* = a* E4 (C a ,T R , f) = C a P T R G E { f) 


(3.8) 

(3.9) 


o l i = a* E4 (C ts , T r . f) = Cfi T r G e ( f) 
and the incoming mean scattering coefficients are given by 
°e* = ** E4 (C a ’ 8 ,T R , f) = Cp s T R G E ( f) 
a\* = a* Fi {CP' s ,T R , f) = c£’ s T r G f ( f) 

&G* = V*Gi(C^ S , Tr, f, SI) = cy Tr Gg( f, SI). 

The integrated RTE, 

-d t J+Sl-VJ= ^aT i -a%*J-a i yj+^-(ca 8 E *ER+a 8 F *F R -St+ca 8 G * V P l R HViV) 

(3.13) 


(3.10) 

(3.11) 

(3.12) 




becomes using the previous definitions : 

-dt J + SI ■ VJ = ^aT 5 - C a P T R G E {f)J - C^T R G E (f)J 

C 47T 

+ j^(cC?,G E (f)E H + C%G F {f) F r • O + cCpG G (f, SI) ^ pyWiV). 
n i , j = 1..3 

(3.14) 

If isotropy is assumed for the functions G EiF}G , the following simpler form is obtained: 
-d t J + Sl-VJ = aT 5 - C a pT R J - C^TrJ 

C 4 7T 

Tr 
4n 


+ -£{cC$Er + C$,F r -SI + cC% £ P l R HViV) 

ij .'J..3 


(3.15) 
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with Cp = 5/4 Cp. In Ripoll & Wray 2004a, it was found that the function G E played an 
important role. Eq. 3.15 should hence not be used for radiating flows where the anisotropy 
||/|| is larger than 0.3. 


4. Macroscopic radiation models for absorbing and scattering media 

We show in this section that the mean coefficients can be easily included in macroscopic 
moment models. We define the three first moments with respect to direction as 

E%{t,r,v) = ^ {I(t, r,Sl,v)) Q (4.1) 

Fg(t, r, i/) = <nj(t,r, «,!/))„ (4.2) 

p "(*> r , 1 ') = i(ft<g>ni(t, r,«,v)) n (4.3) 

The first moment equation is obtained by integrating Eq. (2.1) with respect to f2. We 
obtain 

d t E% + V • Fg = 4nca a {is)B - co a {v)E% (4.4) 

where the incoming and isotropic scattering have canceled using the normalization prop- 
erty (2.4). 

Multiplying (2.1) by and integrating with respect to it, using the phase function 
definition (2.3), we obtain 


Wg + cV • Pg = -(a» + o is {v){ 1 - X (i/)) - o s {v)^- )Fg. (4.5) 

The following three equalities have been used: 

J dSl = 0, J fliflj dn = 47r/3^ j; J dFl = 0. (4.6) 

Using definitions (2.18), (2.19), and (2.20), the integration of (4.4) over frequency gives 

d t E R + V • F R =ca ( a' P T 4 - a a E E R ) (4.7) 

where models a R and a% ~ a%* were given previously in (2.24). 

Integrating (4.5) over frequency leads to the second moment equation 


-dtF R + cV • P R = -{of + ag - \o P f )F r , (4.8) 

C O 

for which the models for the mean coefficients are 

N 

of~of*=Y,^i(C?;T R: f) (4.9) 

i—1 

N 

4 - 4* = E ^i{Ct s ,T R , f) (4.10) 

i = 1 

4~4* = E £7 ^-( (7 f’ s ’ T «’ f ) ( 4 - n ) 

i = 1 


where OFi is defined in (2.39). 

Closure of the macroscopic model (4.7)-(4.8) is achieved by modeling the radiative 
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pressure in (4.8). In many different closures the pressure is written as Pr = DR(f)£jj, 
where D R (f) is the Eddington tensor (Ripoll & Wray 2004b). 

Combining the steady forms of (4.8) and (4.7) to eliminate the flux Fp leads to the 
general Milne-Eddington equations 

“ V ' (7T ! T g V ' (PrE r )) = UpaT A - a%E R . (4.12) 

(a a F + (Jp - P p) 

Or, similarly, by eliminating the energy E R , one obtains 

-V ■ (5« v ■ F«) + (4 + 4 s - \a§.) F R = -cV • (^V R aT 4 ) . (4.13) 

a E 3 a% 

We have two main comments on the model (4.7)-(4.8); they also apply to the formu- 
lation (4.12)-(4.13). First, the normalization property (2.4) eliminates the pressure term 
coming from the scattering in the first moment equation (4.7). Second, the contribution 
of the incoming scattering integrated over direction only enters through the first order 
and delta- function parts of the scattering function in (4.13), since the zero and second 
order terms vanish. 


5. On the use of these models and their numerical costs 

The models presented here will be useful when the solution for frequency dependent 
intensities cannot be done due to its cost. This is usually the case for coupled problems. 
We now give a discussion on reducing the computational cost of the models presented 
here. 

First, when using mean coefficients, it is possible to reduce the computational cost by 
noticing that when the radiation is isotropic and close to equilibrium, the coefficients 
o e,g are equal to the Planck mean and ctj? to (i + l/4)crp. There is hence no need to 
evaluate the complex expressions for oe.f.g-, and they should be simply replaced by their 
limits. More generally, these limits can be extended to ||/|| < 0.1 when T ~ T R . In the 
case where the radiation is isotropic but T ^T R , the limit that should be taken is also 
the Planck mean but evaluated at T R instead of T f. It should also be noticed that the 
form of the mean coefficients proposed here, in terms of A and B, should be retained. As 
a matter of fact, this choice allows checking the different limits and avoids introducing 
singularities. For instance, (2.30) can directly be replaced by their limits, respectivelly 
A = 0 and B = 1/T, for ||/|| < 0.1. 

The solution of the RTE with mean coefficients requires iteration since the mean 
coefficients are nonlinear functions of I. The required number of iterations has been found 
to be small in many simple cases (Ripoll & Wray 2004a), but one could be tempted to 
reduce it further. In the case of using such an equation for coupled problems, the previous 
time step provides excellent starting values for the iterative solution so that the required 
number of iterations should be lower, perhaps only one or two. Another alternative in 
this case could be not to iterate the RTE at all, i.e., to lag it in time, assuming that a 
small difference between the opacity and the intensity due to their non-synchronization 
in time with the fluid motion will not strongly affect the solution. This is likely when a 
global convergence process of the hydrodynamics and radiation to a steady state leads 
to synchronizing all the variables at the end. For the very first iteration, the Planck limit 
seems to be a good initial condition. It should also be noticed here that the moment model 


t as was shown in (3.15) in the linear case 
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(4.7)-(4.8) can similarly be solved by using the mean coefficients from the previous time 
step of the radiation loop or of the hydrodynamics loop in coupled problems. 


6. Conclusion 

The objective of this work was to propose a simple model accounting for radiation 
in complex emitting, absorbing, and scattering media. To do that, models for mean 
absorption and mean isotropic and incoming scattering coefficients have been proposed 
in the case where the various spectral coefficients can be written in terms of polynomial 
functions. Some of these models were previously derived and validated (Ripoll et al. 
2001, Ripoll & Wray 2004a) for non-scattering media; they have been extended here to 
the general case. An integrated RTE which uses these coefficients has been derived where 
the integral scattering term has been absorbed into the modeling. Such a form of the RTE 
is much less costly to solve than the RTE in its non-modeled form. Macroscopic moment 
radiation models, written in their hyperbolic or diffusive forms, have also been derived 
using these coefficients. The particular case where the spectral coefficients are linear in 
frequency has been treated. This case is particularly important for soot and hence for 
combustion applications. It has also been explained how such models can be used at 
a lower cost by reducing the number of iterations needed. We believe the formulations 
proposed here could be used in many complex or coupled problems, like flows radiating 
in dusty media, where, for instance, isotropic and non-isotropic scattering are usually 
disregarded and neglected due to their computational cost. 
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